Accepted for publication in the Astrophysical Journal Letters 
2002 August 8 



On the Low and High Frequency Correlation in Quasi-Periodic 
Oscillations Among White Dwarfs, Neutron Star and Black Hole 

Binaries 

O . ., , 

O ' Lev Titarchuk 1,2 and Kent Wood 2 

3 ■ ABSTRACT 
< 

ON. 

We interpret the correlation over five orders of magnitude between high fre- 
quency Vhigh and low frequency Ui ow in a quasi-periodic oscillations (QPO) found 
by Psaltis, Belloni & van der Klis (1999) for black hole (BH), neutron star (NS) 
! systems and then extended by Mauche (2002) to white dwarf (WD) binaries. The 

observed correlation strongly constrains theoretical models and provides clues to 
C\) '. understanding the nature of the QPO phenomena at large. We argue that the 
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observed correlation is a natural consequence of the Keplerian disk flow adjust- 
ment to the innermost sub-Keplerian boundary conditions near the central object 
which ultimately leads to the formation of the sub-Keplerian transition layer (TL) 
between the adjustment radius and the innermost boundary (the star surface for 
NS and WD and the horizon for BH). In the frameworks of the TL model Phigh is 
related to the Keplerian frequency at the outer (adjustment) radius z/ K and vi ow 
is related to the magnetoacoustic oscillation frequency vma- Using a relation 
between vma the magnetic and gas pressure and the density and the hydrostatic 
equilibrium condition in the disk we infer a linear correlation between z/r and 
Vma, Identification of Vhighi Viow with z/r, vma respectively, lead us to the deter- 
mination of H/r out = 1.5 x 10~ 2 and (5 = 0.1 (where H is the half-width of the 
disk and (3 is a ratio of magnetic pressure to the gas pressure). We estimate the 
magnetic field strength near the TL outer radius for BHs NSs and WDs. The fact 
that the observed high-low frequency correlation over five orders of magnitude is 
valid for BHs, NSs, and down to WDs strongly rules out relativistic models for 
QPO phenomena. We come to the conclusion that the QPOs observations indi- 
cate the adjustment of the geometrically thin disk to sub-Keplerian motion near 
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the central object. This effect is a common feature for a wide class of systems, 
starting from white dwarf binaries up to black hole binaries. 

Subject headings: accretion, accretion disks — stars: neutron — (stars:) white 
dwarfs — stars: oscillations (including pulsations) — X-ray: binaries — (magnetoh- 
ydrodynamics:) MHD 

1. Introduction 

Quasi-periodic oscillations (QPOs) in the bright white dwarf (WD), neutron star (NS) 
and black hole (BH) binaries provide invaluable information on the accretion dynamics in 
the innermost parts of these systems. The question is what kind of objective information 
we can extract from observations using the fundamental principles of fluid physics, radiative 
transfer theory and oscillatory processes? In this sense the discovery of low 20-50 Hz QPOs 
in the luminous NS binaries (van der Klis et al. 1985), and discoveries of NS kilohertz 
QPOs by Strohmayer et al. (1996) and of BH hectohertz by Morgan, Remillard & Greiner 
(1997) opened a new era in the study of the dynamics near compact objects. Psaltis, Belloni 
& van der Klis (1999), hereafter PBK, demonstrated that these NS and BH low and high 
frequencies follow a tight correlation. Their result suggests that the same types of variability 
may occur in both neutron star and black systems over 3 orders of magnitude in frequency. 
These two features are "horizontal branch oscillation" (HBO) along with "low frequency 
noise Lorentzian", ui ow and the lower kHz QPO Uhigh respectively. Belloni, Psaltis & van 
der Klis (2002), hereafter BPK, have updated PBK's correlation adding more points from 
Nowak (2000), Boirin et al. (2000), Homan et al. (2002), Di Salvo et al. (2001) and Nowak 
et al. (2002) data. PBK suggest that the low and high frequencies correlate in a way that 
seems to depend only weakly on the properties of the sources, such the mass, magnetic field, 
or possibly the presence of a hard surface in compact object. At large scales (over 3 orders of 
magnitude) the PBK-BPK correlation for NS and BH systems is fitted by a linear function 
v iow = 0.081///^ and it also has a fine structure at smaller frequency scales (see Fig.l). 

Recently Mauche (2002) has reported low and high frequency QPOs in the dwarf nova 
SS Cygni and VW Hyi (see Wouldt & Warner 2002 for for the VW Hyi observations), and 
has called attention to the fact that these QPOs do extend the correlation of PBK downward 
in frequency by more than two orders of magnitude (see Fig.l). 

In fact, Titarchuk, Lapidus & Muslimov (1998), hereafter TLM98 were first to put forth 
the possibility of dynamical adjustments of a Keplerian disk to the innermost sub-Keplerian 
boundary conditions to explain most observed QPOs in bright low mass X-ray binaries 
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(LMXBs). TLM98 concluded that an isothermal sub-Keplerian transition layer between the 
NS surface and its last Keplerian orbit forms as a result of this adjustment. The TLM98 
model is a general treatment applicable to both NS and black hole systems. The primary 
problem in both NS and BH systems is understanding how the flow changes from pure 
Keplerian to the sub-Keplerian as the radius decreases to small values. TLM98 suggested 
that the discontinuities and abrupt transitions in their solution result from derivatives of 
quantities such as angular velocities (weak shocks). Titarchuk & Osherovich (1999) and 
Titarchuk et al. (2001) interpreted the low noise frequencies as viscous oscillations and 
expected them to provide reliable estimates of radial velocities that would be very close to 
the magnetoacoustic (MA) velocities. The formulation of the MA problem and derivation of 
MA frequencies are presented by Titarchuk, Bradshaw & Wood (2001), hereafter TBW. They 
demonstrate that the observed correlation of the low frequency on the low kHz frequency 
is perfectly described by the dependence of the inferred MA frequency on the Keplerian 
frequency. 

In this Letter, we explain a general correlation between low frequency and high frequency 
oscillations in the TL. In §2, we refer to the QPO data obtained from NSs and BHs by PBK 
and BPK, from the dwarf nova SS Cygni by Mauche (2002) and from VW Hyi by Woudt 
& Warner (2002). The formulation of the problem of low frequency oscillations and the 
relationship with MA oscillations in the transition layer are described in §3. The derivation 
of the low-high frequency correlation is also present in §3. Our summary and conclusions 
follow in §4. 



2. Data Sources 

In Figure 1 we present data points showing low-high frequency dependence from BPK 
(see also PBK), Mauche (2002) and Woudt & Warner (2002) along with a fitted line v low = 
0.081///^. The plot combines data points for BH, NS and white dwarf (WD) systems. A 
reader can find the details of these observations and the QPO data analysis in BPK, PBK and 
Mauche (2002). Each data point is strickly observational, representing two QPO frequencies 
that are detected at the same time in some particular source. 



-4- 



3. High-Low Frequency Correlation 

3.1. The main features of the transition layer model 

TLM98 define the transition layer as a region confined between the the inner sub- 
Keplerian disk boundary and the first Keplerian orbit (for the TL geometry, see Fig. 1 in 
TLM98 and Fig. 1 in Titarchuk, Osherovich & Kuznetsov 1999). The radial motion in the 
disk is controlled by friction and the angular momentum exchange between adjacent layers, 
resulting in the loss of the initial angular momentum by accreting matter. The corresponding 
radial transport of the angular momentum in a disk is described by the following equation 
(see e.g. Shakura & Sunyaev 1973): 

M^M 2 ) = 2n^L(W rv R 2 ), (1) 

where M is the accretion rate in the disk and W rip is the component of a viscous stress tensor 
that is related to the gradient of the rotational frequency uj = 2ttu, namely, 

W rv = -2vHR^, (2) 

where H is a half-thickness of a disk and rj is turbulent viscosity. The nondimensional 
parameter that is essential for equation (1) is the Reynolds number for the accretion flow, 

7 = M /A-KTjH = RV r /D, (3) 

where V r is a characteristic velocity, and D is the diffusion coefficient. D can be defined as 
D = V t l t /3 using the turbulent velocity and the related turbulent scale, respectively or as 
D = Dm = c 2 /a for the magnetic case where a is the conductivity (e.g. see details of the 
D-definition in Lang 1998). Equations oo = ujq at R = ri n (the inner disk radius), uj = wk 
at R = r out (the radius where the transition layer adjusts to the Keplerian motion), and 
du/dr = duj K /dr at R = r out were assumed by TLM98 to be these boundary conditions. 
Thus, the profile uj(R) and the outer radius of the transition layer are uniquely determined 
by the boundary conditions and the angular momentum equation (1, 3) for a given value 
of Reynolds number 7. TLM98 present the analytical solution of this TL problem (see also 
T099 , Eqs. 8-9). For the observed range of v k q L = 400 — 1000 Hz one can find using the 
TLM98 solution that 7 varies from one to 5 (see also Fig. 3 in TLM98) . 



3.2. Magnetoacoustic oscillations in the disk transition layer 



The frequency of the QPO associated with MA oscillations and the correlation of the 
MA frequency with the Keplerian frequency u K were derived by TBW. The MA frequency is 
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derived as the eigenfrequency of the boundary-value problem resulting from a MHD treat- 
ment of the interaction of the disk with the magnetic field. The problem is solved for two 
limiting boundary conditions encompassing realistic possibilities. The solution yields a veloc- 
ity identified as a mixture of the sound speed and the Alfven velocity. The TBW treatment 
does not specify how the eigenfrequency is excited or damped. However, it makes clear that 
the QPO is a readily stimulated resonant frequency. TBW construct an approximate formula 
for the MA frequency i>ma using the asymptotic forms for the two extreme (acoustic and 
magnetic) cases (see TBW, Eqs. 12-13) 

Vma « l(b s /n) 2 v 2 s + (&mA)^m] V2 , (4) 

where bu = b s = n for the stiff boundary conditions, bu and b s are determined by Eqs. 
(16-17) in TBW for the free boundary conditions. They are approximately bu = n and 
b s = 1. The frequencies v a , v§ — vm (cx — 6) and u — v s (a — 0) are determined by the 
equation 

v a = [(a + 2)/A]AW/(r<£W-r™<% (5) 

where A = A s = s 2 and s 2 is the square of the sound velocity for the pure acoustic case, and 
a = 6 and A = A m = B 2 o ^ t r" ut /4:7ip for the dipole magnetic case (p is the TL mean density). 
For the case of interest r out /r in > 2 we approximate 

v s = s/2(r out - r in ) ~ s/2r out (6) 

and 

VM = (WY'HrT-t(r,Jr m ,n ~ ™*~)/*~ < 7 > 
where Vx{r ou t = B(r out )(4:TTp)^ 1 ^ 2 is the Alfven velocity at r out . 

Now we express the MA oscillation frequencies through the gas pressure P g = s 2 p and 
the magnetic pressure P M = V%p/2 using equations (4,6,7) 



v 



.2 

MA 2 

P r out 



(fP 9 + PM), (8) 



where / = 1/32, l/(327r 2 ) for the stiff and free boundary conditions respectively. 



3.3. Low-high frequency correlation 

If j3 — Pm/P 9 is a ratio of the magnetic pressure to the gas pressure we can rewrite 
equation (8) as follows 

^J3— (9 > 
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The hydrostatic equilibrium in the disk is presented by the equation 

P/p = (P g + Pm)/ P = (2nu K ) 2 r 2 out (H/r out ) 2 . (10) 
Using Eq. (9) we can transform this relation into 

= {2^ K frUH/r out ) 2 . (11) 



v 1 + (3) l J \i : \l'i„, _ _ ■) . 2 



which leads us ultimately to the linear relation between uma and vk 

vma = Cma^k- (12) 

where C M a = 2 1//2 47r[(/ + f3)/{\ + / 9)] 1 / 2 (i/ /r out ) is a proportionality coefficient which will 
be "universal" to the extent that f3 and H/r ou t remain about the same from one source to 
the next. To proceed further we estimate the magnitude of Cma separately for NS, BH and 
WD systems. Hydrostatic equilibrium (Eq.10) gives 

H/r out (P g / p) 1/2 

S/V K \ r=rout , [id) 



(1 + /3)V2 v K (r out ) 

where Vk = (GM/r) 1 / 2 is the Keplerian velocity. Approximate constancy within a source 
class of the coefficient Cma oc H/r out = s/Vk can perhaps be understood as the thermaliza- 
tion of an excess kinetic energy according to the virial theorem. If the disk is in equilibrium 
with nonradial inward motion then all of the orbital velocity at r out would be needed for virial 
equilibrium at that radius. The small amount of inward motion (see Shakura & Sunyaev 
1973) creates a disposable energy budget given U excess = m p V 2 /2. Thus one can conclude 
that (Uexcess/Ugrav)) 1 / 2 sets the maximum value of H/r out . After that, the argument is that 
it grows to that maximum because there is little else one can do with the excess except 
thermalize. The coefficient Cma applicable within a class is estimated from canonical or 
representative values for the Keplerian frequency (essentially the lower kHz QPO frequency 
< u kQL >), the plasma temperature T, all evaluated at the transition layer outer boundary 
r out and the scaled central object mass, m = M/M Q (in the solar units). We obtain 

Cma = 0.08 x kf^ 2 m w 1/3 < v kQL ^ >~ x ^ (/?/0.1) 1/2 , (14) 

where A;T 25 = kT /2b keV, m w = m/10 and < UkQL,i5o >=< VkQL > /150 Hz. 

Laurent and Titarchuk (1999) demonstrated that the temperature of the Compton 
cooled corona kT is about 25 keV and 12 keV for BH and NS respectively. This assumes 
that the energy released in the corona (supposedly the TL) is comparable with the external 
illumination flux, coming from the disk and from the star surface (in the NS case). Mauche 
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(2002) and Woudt & Warner (2002) detected the QPO frequencies from white dwarfs (WD), 
SS Cyg and VW Hyi during outbursts when the temperature of the transition layer should 
be very close to the effective temperature of the star surface, i.e. kT = 0.05 — 0.1 keV. With 
assumptions kf = 25, 12, 0.05 keV, m = 10, 1.4, 1, < v kQL >= 150, 800, 0.1 Hz for 
BH, NS and WD respectively and f3 = 0.1 we obtain Cma values that are very close to the 
observed one, C MA ~ 0.08. They are 0.08, 0.061 and 0.089 for BH, NS and WD respectively. 
Formula (13) also gives us estimate of the ratio H/r out using kT, m and u k Q L 

H/r out = 0.015 x 4 /2 m~ 1/3 < v kQL , 150 >" 1/3 . (15) 

Thus H/r out « 0.015, 0.011, 0.017 for BH, NS and WD respectively. 

3.4. B-field estimate 

Now we can proceed with the magnetic field strength determination in the TL using the 
obtained value of (3 = 0.1. Because Pm = (3P g , P g = (kT/m p )p and pm = B 2 /8tt we find 
that 

B = 1.14 x lO^fm^V/ 2 < v kQL , 150 > x / 3 G, (16) 

where r = {p/m p )a^r out is the Tompson optical depth of the TL. Thus 5 ^ 1.1 x 10 5 , 2 x 
10 5 , 700 G in the TL of BH, NS and WD systems respectively. 

4. Summary and Conclusions 

We have presented a model which treats the radial oscillations in the transition layer 
surrounding a black hole, neutron star and white dwarf. 

(1) We found that the observed high- low frequency correlation is a natural consequence 
of an adjustment of the Keplerian disk flow to the innermost sub-Keplerian boundary con- 
ditions near the central object. This ultimately leads to the formation of the sub-Keplerian 
transition layer (TL) between the adjustment radius and the innermost boundary (the star 
surface for NS and WD and the horizon for BH). (2) In the framework of the TL model Vhigh 
is related to the Keplerian frequency at the outer (adjustment) radius vk and v\ ow is related 
to the magnetoacoustic oscillation frequency uma- Using a relation between uma the mag- 
netic and gas pressure and the density and the hydrostatic equilibrium condition in the disk 
we infer a linear correlation between z/ K and u M a (see formula 12). (3) Identification of v high 
and vi ow with uk and uma respectively, lead us to the determination of H/r out ~ 1.5 x 10~ 2 . 
(4) We present strong arguments that the ratio (3 = Pm/P 9 ~ 0.1 in all binaries studied. (5) 
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We estimate the magnetic field strength near the TL outer radius for BHs, NSs and WDs. 
(6) The fact that the observed correlation of high and low frequencies holds over five orders 
of magnitude shows that the QPO frequencies cannot be set by general relativistic effects, 
because they are essentially non-relativistic in the white dwarf case. 

We acknowledge Chris Mauche and Tomaso Belloni for kindly supplying us data shown 
in Figure 1 and Chris Shrader, Paul Ray and particularly, Chris Mauche for fruitful discus- 
sions that motivated this paper. 
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Fig. 1. — Correlation between frequencies of QPO and noise components white dwarfs [dia- 
monds for SS Cyg (Mauche 2002); squares for VW Hyi (Woudt & Warner 2002)] , neutron 
star (open circles) and, black hole candidate (filled circles) sources. Neutron star and black 
hole data are from Belloni, Psaltis & van der Klis (2002). The dashed line represents the 
best-fit of the observed correlation. This plot also appears in Mauche (2002). 



